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Abstract:  The  properties  of  the  flow  near  the  trailing  edge 
of  a  body  of  revolution  with  fins  has  been 
investigated  experimentally  and  compared  with 
predicted  values  baaed  on  existing  thewies,,  The 
theoretical  predictions  were  based  on  the  solution 
to  the  axlsymmetrlc  boundary-layer  equations  of 
Truckenbrodt  and  the  flat-plate  solution  of 
Von  girnuin  and  Pohlhausen.. 

Experimental  work  included  measurements  with  a 
hot-wire  anemometer  to  determine  same  statistical 
properties  of  the  randomly  fluctuating  velocity 
field 0  These  properties  include  correlation 
functions  and  their  related  power-spectral  densities 
and  steady  mean  values  and  root- me an- square  values 
of  the  velocity  f luctuatiocsa 

Further  properties  of  the  flow  field  were  established 
with  a  pitot-static  tube  and  a  vorticity  meter  that 
measured  streanmise  vortlcity  arising  from  the  fin- 
and-body  configuration..  Measurements  were  also  made 
for  the  bf'dy  with  a  propeller  located  at  0  97  per 
cent  of  the  body  length  from  the  nose  The  results 
obtained  establish  that  the  aft  fins  are  the  major 
cause  of  the  radial  and  circumferential  fluctuation 
of  the  mean  vel  x^ity  which  give  rise  to  unsteady 
forces  in  the  propulsor, 
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CHAPTER  I 


INTRODUCTION 

A-  Discussion  of  the  General  Field  of  Study 

The  flow  behind  a  body  of  revolution  with  fins  is  of 
practical  importance  in  aeronautical  and  marine  engineering- 
Although  the  flow  of  real  fluids  over  bodies  of  revolution  and  over 
airfoils  are  well  understood,  the  fluid  motion  over  a  body  of 
revolution  equipped  with  fins  cannot  be  predicted.  In  particular, 
the  effects  of  viscosity  on  the  flow  at  the  intersection  of  a  fin 
and  a  cylinder  (or  even  a  flat  plate)  still  awaits  a  theoretical 
explanation.  It  should  be  noted  that  this  effect  largely  determines 
the  characteristics  of  the  flow  field  in  which,  for  instance,  the 
propulsor  of  an  underwater  vehicle  operates*  In  such  an  application, 
the  designer  requires  an  accurate  knowledge  of  the  inflow  velocity 
field  to  be  able  to  establish  the  geometry  of  a  propulsor  capable  of 
developing  the  necessary  thrust  with  a  minimum  level  of  unsteady 
forces- 

The  solution  to  the  above-mentioned  problem  can  be  obtained 
from  the  solution  to  its  separate  parts: 

a)  In  the  first  stage,  the  properties  of  the  flow  near  the 
trailing  edge  of  a  body  of  revolution  with  fins  must  be 
established  experimentally  and  compared  with  predicted 
values  based  on  existing  theories-  This  information  will 


£ 


reveal  the  true  nature  of  the  flow  and  the  limitations 
of  the  theoretical  approaches. 

b)  In  the  second  stage,  detailed  experimental  and 
theoretical  studies  of  special  flow  problems  must  be  made 
in  order  to  explain  the  discrepancies  observed*  As  an 
example,  the  flow  of  a  real  fluid  near  the  Intersection 
of  a  strut  and  a  flat  plate  must  be  studied. 

c)  In  the  third  stage,  attempts  must  be  made  to  control  the 
nature  of  the  wake  by  practical  or  technically  simple 
means  that  can  he  used  on  existing  vehicles- 

B.  Statement  of  the  Problem 

The  primary  objectives  of  this  thesis  are  those  stated  in 
paragraph  a).  In  particular  it  was  proposed  to  measure  various 
properties  of  the  flow  near  the  trailing  edges  of  a  slender  body  of 
revolution  with  fins  and  to  compare  these  measurements  with 
theoretically  predicted  values-  In  addition,  the  effects  of  a 
propeller  on  these  observations  was  assessed-  More  precisely,  the 
experimental  work  included  the  following  measurements; 

1-  Circumferential  and  radial  variations  of  mean  velocity 
and  the  associated  harmonic  content 

2.  Boundary- layer  profiles 

3-  Circumferential  and  radial  variations  of  the  random 
velocity  fluctuations 

4o  Auto-  and  cross-correlation  functions  of  the  random 
velocity  fluctuations 

5-  Spectral  and  cross-power  spectral  densities  of  the  random 


velocity  fluctuations 

6-  The  detection  and  measurement  of  stream-wise  vortices 
originating  at  the  fin-and-body  intersections 

7-  A  preliminary  assessment  of  the  effects  of  blowing  at  the 
fin-and-body  intersections* 

The  theoretical  work  consisted  of  predictions  of  the 
following  quantities; 

1*  The  velocity  field  for  a  streamlined  body  of  revolution 
without  fins  in  a  real  fluid 

2.  The  velocity  field  behind  the  fins 

J.  The  interaction  between  body  and  fins. 
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CHAPTER  II 

SUMMARY  OF  THE  INVESTIGATION 

The  flow  near  the  trailing  edge  of  a  slender  body  of 
revolution  with  fins  has  been  investigated  by  means  of  three 
different  instruments:  Pitot  tubes,  hot-wire  anemometers,  and  a 
vorticity  meter.  The  measurements  were  made  in  an  axial  plane  at  a 
distance  from  the  leading  edge  of  the  model  equivalent  to  95  per 
cent  of  the  body  length,  i.e»,  in  a  region  normally  occupied  by  the 
propulsor  of  a  typical  underwater  vehicle.  The  mean -velocity 
profiles  as  well  as  some  statistical  properties  of  the  random 
velocity  fluctuations  were  established-  In  addition,  the  streamwise 
vortices  formed  at  the  body-and-fin  intersections  were  investigated- 
The  experiments  were  performed  in  a  wind  tunnel  at  a  Reynolds  number 
of  5.5  x  10  ,  based  on  body  length- 

The  steady  (mean)  flow  observed  near  the  body  showed  little 
circumferential  and  radial  variations  when  the  body,  supported  by  a 
single  airfoil  strut,  had  no  aft  fins.  Figure  1  shows  a  typical 
part  of  the  circumferential  distribution  of  mean  velocity,  and 
Fig-  2  shows  corresponding  harmonic  components.  When  the  four  aft 
fins  were  attached  to  the  body,  just  upstream  of  the  plane  of  the 
measurements,  the  disturbances  of  the  mean  flow  became  more 
pronounced,  as  shown  in  Figs-  1  and  2- 

It  was  clear  from  a  comparison  of  calculated  and  measured 
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Fig.  2  Harmonic  Content  for  Various  Body  Configurations 
at  r/D  =  0.10 
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mean-velocity  profiles  (Fig-  3)  that  the  aft  fins  affect  the  flow 
in  two  ways-  At  relatively  large  distances  from  the  body  axis*  such 
as  at  r/D  -  0-35*  sharp  depressions  in  the  velocity  profile  are 
produced  that  agree  with  theoretical  predictions-  Near  the 
body-and-fin  intersection,  however,  vortices  were  formed  as  shown 
in  Fig-  4=  These  vortices  entrain  fresh  fluid  into  the  immediate 
neighborhood  of  the  fins  and  displace  sluggish  fluid  from  the 
boundary  layer-  As  a  result,  large  peaks  and  valleys  were  created 
in  the  velocity  profiles  (refer  to  curve  for  r/D  ®  0-10  in  Fig-  3)° 
Thus,  Figs-  1  to  k  establish  conclusively  that  the  aft  fins  are  the 
major  cause  of  the  radial  and  circumferential  fluctuations  of  the 
mean  velocity,  which  often  produce  strong  oscillatory  forces  in  the 
propulsors- 

Before  concluding  this  investigation  an  attempt  was  made  to 
"de-spin"  the  corner  vortices  by  blowing  air  through  a  pair  of  small 
holes  located  near  the  leading  edge  of  one  of  the  fins  (Fig-  5)- 
The  size  and  location  of  the  holes  were  arbitrarily  chosen; 
therefore,  the  resulting  conditions  are  not  to  be  considered 
optimum-  Nevertheless,  the  vorticity  probe  (Fig*  6)  indicated  a 
marked  decrease  in  the  streamwise  rotation.  The  amplitudes  of  the 
mean  velocity  variations  were  also  reduced,  as  shown  in  Fig-  7- 
This  method,  therefore,  appears  to  hold  promise  for  the  eventual 
control  of  the  velocity  field  in  which  propulsors  operate,  and 
further  investigations- -both  theoretical  and  e,xperimental~-will  be 
made- 

Measurements  made  with  two  constant- temperature  hot-wire 


on  Each  Side  of  Fin-Body  Intersection 
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AIR  BLOWN  NORMAL  TO  FREE-STREAM 
AT  ANGLE  4>  =  45  deg  BETWEEN  STRUT 
AND  PLATE 


Fig.  5  Location  of  Holes  for  Blowing  Air  to  Reduce  Intensity 
of  Vortex 
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Fig.  6  Vorticity  Probe 


anemometers  established  a  number  of  statistical  properties  of  the 
turbulent  velocity  fluctuations*  Figure  8  shows  a  typical 
circumferential  distribution  of  the  root-mean- square  values*  During 
the  tests,  these  values  varied  from  a  minimum  of  3  per  cent  to  a 
maximum  of  10  per  cent.  The  power  spectra  corresponding  to  various 
test  conditions  are  shown  in  Fig.  9°  The  peak  energy  occurs  at  a 
nondimensional  frequency  where  the  characteristic  length  is  the 
maximum  diameter  of  the  body.  Finally,  Figs.  10  to  13  show  typical 
autocorrelation  and  cross- correlation  functions,  as  well  as 
cross-power  spectral  densities*  These  were  calculated  from  magnetic 
tape  recordings  made  during  the  tests  with  the  aid  of  special 
analog- to-digital  converters  and  high-speed  digital  computers.  Most 
of  the  measurements  were  made  both  with  propeller  operating  at  its 
design  advance  ratio  and  with  no  propeller. 

The  effects  of  a  propeller  are  most  notable  in  .Fig.  8°  The 
propeller  creates  a  favorable  pressure  gradient  at  a  critical 
position  on  the  body  thereby  reducing  the  mean,  velocity  fluctuations 
as  well  as  the  intensity  of  the  turbulent  motions  of  the  flow. 

Figure  2,  which  shows  a  reduction  in  the  harmonic  content, 
illustrates  some  of  the  beneficial  influences  produced  by  the 
propeller |  however,  it  is  difficult  to  distinguish  the  influences 
produced  by  the  induction  effects  from  those  resulting  from 


viscosity 


3A  3avrffis-NV3i/\i-±ooa 


Fig.  8  Typical  Root-Mean-Square  Values  of  the  Random-Velocity 
Fluctuations  for  Body  with  Fins 
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Fig.  11  Typical  Normalized  Cross-Correlation  Function 


90000 


Fig.  12  Real  Part  of  Cross-Power  Spectral  Density  Associated 
with  Fig-  11 


0.0006 


Associated  with  Fig.  11 


CHAPTER  III 


EXPERBffiMTAL  STUDY 

A-  Discussion  of  Test  Procedure  and  Alms 

The  experimental  measurements  were  made  in  the  wake  of  a 
streamlined  body  of  revolution  with  fins.  The  body  configuration 
approximated  an  8:1  ellipsoid  of  revolution  as  shown  in  Fig.  14. 

The  main  objective  of  the  investigation  was  to  obtain  the 
flow  characteristics  at  the  region  of  inflow  to  the  propeller.  The 
measurements  were  made  at  various  axial  and  radial  locations  as  well 
as  circumferential  positions.  Steady- state  components  were  obtained 
from  which  mean-velocity  distributions  were  obtained.  Measurements 
were  also  made  of  the  root-mean- square  (rms)  values  of  the  randomly 
fluctuating  velocity  components.  This  resulted  in  the  determination 
of  the  turbulence  level  of  the  flow  field.  Auto-  and  cross- 
correlation  functions  were  calculated  on  the  basis  of  measurements 
taken  at  different  points  in  the  flow  field-  Inherent  in  the 
reduction  of  the  correlation  data  are  several  mathematical 
techniques t  these  consisted  of  performing  Fourier  Transforms  to 
obtain  normalized  power  spectra  and  Fourier  Analyses  to  obtain 
harmonic  content.  The  normalized  power  spectra  were  also  obtained 
directly  from  the  measurements  by  recording  the  signals  from  the 
measuring  device  on  magnetic  tape  and  playing  them  back  through  a 
wave  analyzer  and  integrator. 
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The  tests  were  conducted  in  the  48-inch,  diameter*  16-foot 
long  octagonal  test  section  of  the  wind  tunnel  at  the  Garfield 
Thomas  Water  Tunnel*  Ordnance  Research  Laboratory*  at  The 
Pennsylvania  State  University.  For  this  investigation*  the  air 
speed  in  the  tunnel  was  maintained  at  approximately  100  fps»  A 
detailed  description  of  this  closed-circuit  tunnel  will  be  found  in 
Reference  1. 

To  maintain  a  complete  but  efficient  test  program*  the 
position  variables  were  kept  at  a  minimum-  The  measurements  were 
made  at  a  constant  axial  location,  between  the  control  surface 
trailing  edge  and  the  propeller  leading  edge.  The  radial  positions 
were  varied  outward  from  the  body  while  circumferential  traverses 
were  made  completely  around  the  body. 

To  present  the  results  in  as  general  a  manner  as  possible* 
certain  constants  were  chosen  by  which  the  variable  quantities 
were  nondimensionalized-  Thus,  linear  dimensions  were  non- 
dimensi Dualized  with  respect  to  body  length*  body  diameter,  or 
displacement  thickness.  Velocity  measurements  were  normalized  with 
respect  to  the  free-stream  value  or  some  other  characteristic 
parameter,  such  as  the  mean  velocity  in  the  boundary  layer. 

The  primary  independent  variables  were  the  radial  and 
circumferential  positions  around  the  body.  The  radial  positions 
chosen  were  at  values  for  r/D  of  0.10*  0.17,  0.24,  and  Q«35J 
where  r  is  the  normal  distance  from  the  body  axis  to  the  point  of 
interest  and  D  is  the  maximum  body  diameter-  The  circumferential 
traverses  were  performed  completely  around  the  body  for  each  value 
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of  r/D. 

The  measurements  were  performed  without  the  propeller  and 
then  repeated  after  the  propeller  was  added;  thus,  any  difference 
in  the  measurements  was  due  solely  to  the  presence  of  the  propeller. 

Figure  15  shows  the  radial  positions  of  the  measurements  and 
the  coordinate  system  chosen  to  describe  the  circumferential 
traverses.  All  measurements  were  performed  with  a  constant- 
temperature  hot-wire  anemometer  and  checked  against  a  pitot-static 
probe  that  was  used  as  a  reference-  The  hot-wire-anemometer 
equipment  and  its  use  are  described  in  a  later  section-  The 
repeatability  of  the  data  was  within  3  to  5  per  cent  from  day  to  day 
and  from  week  to  week- 

The  model  is  approximated  by  an  8 si  ellipsoid  of  revolution 
with  fins  (Fig.  14).  The  model  was  supported  in  the  tunnel  by  the 
large  fin  at  the  forward  part  of  the  body  rather  than  sting-mounted 
to  eliminate  any  vibrations  that  would  occur  as  a  result  of 
cantilevered  mounting. 

Referring  to  Fig.  15,  the  lower  strut  was  designated  zero 
angle;  positive  angles  were  those  measured  counterclockwise  when 
viewed  upstream-  The  other  fins  were  located  on  the  aft  part  of  the 
body  and  spaced  equidistantly.  The  profile  of  all  struts,  fore  and 
aft,  was  an  HAGA-OOlJ  airfoil  section. 

To  eliminate  any  effects  produced  by  tunnel-wall  interference, 
a  "liner"  was  used  in  the  test  section.  The  resulting  inner  contour 
in  the  test  section  was  determined  by  a  potential  flow  solution  for 
the  body  and  approximated  a  stream  surface.  This  procedure  is 
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necessary  when  the  ratio  of  body  to  test-section  diameter  becomes 
greater  than  1;6°  (in  this  case  the  ratio  is  about  lsj).  As  a 
result  of  this  modification,  the  flow  field  in  the  tunnel  was  almost 
axi symmetric. 

Two  hot-wire  probes  were  used  for  the  velocity  measurements. 
One  probe  was  placed  at  a  constant  axial  position  along  the  body  and 
at  a  constant  angular  position  between  two  control  surfaces-  This 
probe  was  moved  only  in  a  radial  direction-  The  second  prohe  was 
at  the  same  axial  station,  but  its  position  was  varied  in  both  the 
radial  and  circumferential  directions.  This  prohe,  referred  to  as 
the  "movable"  probe,  was  mounted  on  a  traversing  mechanism  attached 
to  the  external  surface  of  the  tunnel  wall. 

The  moveable  probe  was  mounted  on  a  strut  that,  in  turn,  was 
mounted  on  a  platform.  The  platform  was  mounted  on  bars  that  passed 
completely  through  the  tunnel  such  that  it  could  slide  along  the 
bars  to  various  positions  behind  the  body-  The  strut  could  be 
rotated  by  an  external  control-  Thus,  the  second  probe  could  be 
fixed  at  a  radial  position  and  rotated  to  make  the  required 
circumferential  measurements. 

The  signals  from  the  probes  were  recorded  simultaneously  on 
magnetic  tape.  These  signals  were  later  played  back  and  re-recorded, 
for  a  computer  that  introduced  the  proper  time  lags  so  that 
correlation  functions  could  be  obtained. 

The  testing  procedure  was  the  same  both  with  and  without  the 
propeller  Installed  on  the  model.  During  tests  with  the  propeller 
in  place,  all  measurements  were  made  at  its  design  advance  ratio. 


For  these  tests  the  moveable  probe  -was  modified  so  that  measurements 
could  be  made  at  the  inflow  region.  The  probe  holder  was  extended 
so  that  the  probe  was  in  a  vertical  position  in  front  of  the 
propeller  instead  of  horizontal.  This  new  orientation  of  the  probe 
did  not  affect  the  ability  of  the  wire  to  sense  axial  velocity, 
because  it  will  respond  only  to  fluctuations  normal  to  its  own  axis. 
The  position  of  the  extension  device  was  such  that  it  did  not 
interfere  with  the  operation  of  the  propeller  or  disturb  the  flow 
field. 

B-  Test  Equipment 

The  primary  piece  of  test  equipment  was  a  constant  - 
temperature  hot-wire  anemometer  designed  by  Kovasznay  and 
constructed  at  the  Johns  Hopkins  University. (2,3) 

Brief ly5  the  hot-wire  anemometer  is  used  to  measure  random 
velocity  fluctuations  by  producing  a  voltage  that  is  proportional  to 
velocity.  The  wire  (usually  Q.000l5-in. -dia  tungsten)  is  placed  ia 
a  resistance  bridge  and  bridge  balance  is  maintained  by  closing  the 
feedback  loop  with  a  d-c  amplifier.  The  air  flowing  over  the  wire 
cools  it,  thus  changing  its  resistance  and  unbalancing  the  bridge. 
The  unbalance  produced  by  the  air  flow  is  recorded  as  a  fluctuating 
voltage  about  some  predetermined  mean- 

The  physical  law  that  describes  the  phenomenon  of  heat  loss 
from  the  wire  is  assumed  to  follow  King’s  law: 

H  =AT(ATTT  +C) 

where  H  the  heat  loss  of  the  wire,  AT  the  difference  between  the 
temperature  of  the  unheated  wire  and  the  actual  operating 
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temperature  of  the  -wire,  U  the  air-flow  velocity,  and  A  and  C  are 
constants  that  are  determined  by  calibration.  In  the  equilibrium 
stage  of  operation,  this  heat  loss  is  equal  to  the  electrical  power 
input. 

The  appropriate  nonlinear  circuit  to  invert  King's  law 
consists  of  two  squaring  circuits  and  an  adjustable  bias.  The 
circuit  device,  called  a  "linearizer,"  converts  the  output  of  the 
hot-wire  anemometer  into  a  voltage  that  is  proportional  to  velocity. 
The  amplifier  and  linearizer  are  shown  in  Fig.  16. 

Figure  17  is  a  block  diagram  of  the  amplifier  and  linearizer 
with  the  peripheral  equipment  used  for  recording  data-  Table  I 
lists  this  equipment. 

Appendix  I  outlines  a  method  for  checking  the  response  of 
the  hot-wire  anemometer  and  comparing  this  response  with  that  of  a 
pitot-static  probe. 

C-  Method  for  Obtaining  and  Presenting  Test  Results 

The  testing  procedure  is  described  here  in  detail  to  show  the 
evolution  of  the  work. 

1.  Body  with  fins  and  without  propeller 

a.  Mean-velocity  distribution 

b.  Root-mean- square  velocity  distribution 

c.  Boundary- layer  profiles 

2.  Body  with  fins  and  propeller 

a-  Mean-velocity  distribution 

b.  Root -mean- square  velocity  distribution 


c.  Boundary- layer  profiles 


Fig.  17  Block  Diagram  of  Peripheral  Equipment  Used  with  the 
Hot-Wire  Anemometer 


TABLE  I 
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-2- 5V  D-C  Power  Supply  (two  1°5V  dry  cells  with 
10, 000- ohm  potentiometer) 


d.  Boundary- layer  profiles  at  x/L  ~  O.99 
Unless  otherwise  noted,,  the  velocity  distributions  were  made  at 
radial  positions  of  r/D  equal  to  0.10,  0.17,  0.24,  and  0.J5J  all 
circumferential  measurements  were  made  for  0  between  0  and  36O  deg; 
all  axial  positions  remained  constant  at  x/L  -  0-95j  and  the 
boundary- layer  profiles  were  made  at  x/L  *-  0-95  and  0  =s  45  deg. 

The  output  signal  from  the  hot-wire  anemometer  was  fed  into 
a  d-c  voltmeter,  which  recorded  the  mean  value  of  the  fluctuations. 
This  voltage  was  compared  with  the  voltage  found  for  the  free- stream 
velocity  from  calibration  test  runs.  Thus,  at  the  particular  point 
in  question,  the  ratio  of  U  local  to  Uco  was  obtained. 

The  fluctuating  component  was  recorded  on  a  true  xtms  voltmeter. 
When  this  voltage  was  compared  with  the  free-stream  value,  the 
turbulence  level,  or  rms  value,  was  obtained.  Figures  18  and  19  are 
typical  of  the  graphical,  presentation  made  of  these  results.  The 
complete  graphical  description  of  the  flow  field  for  each  test 
condition  is  given  in  Appendix  III. 

The  boundary- layer  profiles  were  obtained  in  the  same  manner 
as  the  wake  profiles,  except  that  the  angle  of  traverse  was  kept 
constant  and  the  radial  position  was  varied.  A  typical  profile  of 
mean  velocity  vs  distance  from  the  body  is  shown  in  Fig.  20. 

The  boundary- layer  characteristics  at  some  axial  station 
are  based  on  the  displacement  and  momentum  thicknesses  rather  than 
just  the  boundary- layer  thickness. (4)  This  allows  for  more 
meaningful  results  to  be  concluded  from  a  particular  investigation. 

The  displacement  thickness,  defined  as  the  decrease  in  volume 
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flow  resulting  from  the  influence  of  friction,  is  given  by 

v  ■$-&)*>  • 
y -o 

where  the  boundary- layer  profile  is  integrated  from  the  body  surface 
to  the  free- stream- 

The  momentum  thickness,  which  represents  the  loss  of 
momentum  in  the  boundary  layer,  is  given  as 

y>  =  f-k  /i-_i  \  j. 


y=o 
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These  quantities  define  the  shape  factor  H  =  O  /'fi  »  which 
describes  the  boundary-layer  profile  at  any  point  along  the  body-  A 
more  detailed  discussion  regarding  the  actual  calculations,  as  well 
as  the  values  of  S  and  2^,  can  be  found  in  Appendix  II. 

B-  Bata-Reductlon  Techniques 

An  oscilloscope  and  a  Panoramic  Sonic  Analyzer  were  employed 
as  visual  aids  for  detecting  turbulence.  The  analyzer  displays  the 
energy  level  versus  frequency  at  a  particular  instant.  This 
analyzer  d.oes  not  give  a  true  power  spectrum;  it  simply  indicates 
regions  of  high  energy  concentrations  that  should  be  more  fully 
investigated. 

The  normalized  power  spectrum  was  obtained  by  passing  the 
hot-wire  signal,  through  a  wave  analyzer  and  integrating  it  for 
approximately  JO  sec-  Figure  21  is  a  comparison  of  this  power 
spectrum  with  that  obtained  by  performing  a  Fourier  transform  of  the 


Fig.  21  Comparison  of  Normalized  Power  Spectra 


autocorrelation  function  shown  in  Fig.  10.  This  last  method 
represents  a  mathematical  technique  for  obtaining  the  normalized 
power  spectra  and  shows  the  results  of  this  method  to  be  in 
agreement  with  those  obtained  by  direct  measurement. 

This  normalized  curve  is  shown  (dotted)  in  Fig.  21.  lote 
that  the  abscissa  in  this  figure  has  been  normalized  with  respect  to 
the  displacement  thickness  and  free- stream  velocity.  This  has  been 
done  to  make  the  results  compatible  with  the  work  of  other 
investigators • ( 5 ) > ( 6 ) 

The  cross-power  spectral  density  was  obtained  in  a  manner 
analogous  to  the  foregoing  except  that  the  cross-correlation  rather 
than  the  autocorrelation  function  is  used-  Figure  LI.  shows  a 
typical  cross- correlation  function  related  to  the  autocorrelation 
function  of  Fig.  10. 

Figures  12  and  13  are  the  cross-power  spectral  densities  of 
the  real  and  imaginary  parts,  respectively,  of  Fig.  11.  From  each 
of  the  figures  showing  the  power  spectrum,  it  is  seen  that,  as  the 
frequency  S*/^  increases,  the  available  energy  approaches  zero. 
This  is  eacpected  since  only  a  finite  amount  of  energy  is  available 
in  the  flow  field. 

The  cross-power  spectral  densities  associated  with  specific 
points  in  the  flow  -field  are  shown  in  Appendix  III  together  with  the 
corresponding  cross- correlation  functions. 

Of  vital  interest  to  the  propeller  designer  is  the  harmonic 
content  of  the  flow  field.  As  an  illustration  of  the  harmonic 
content,  consider  Fig.  22;  this  shows  a  large  amplitude  at  the 
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fourth  harmonic,  which  is  produced  by  the  fins* 

Figure  23  is  a  similar  analysis  for  the  body  with  only  the 
forward  supporting  strut.  As  expected,  there  is  no  real  predominant 
harmonic,  even  close  to  the  body.  The  differences  between  the  two 
figures  is  due  solely  to  the  effects  of  the  afterbody  struts. 

These  graphs  of  harmonic  content  were  obtained  by  performing 
a  Fourier  Analysis  on  the  wake  velocity  profiles  for  various  values 
of  r/D*  (These  profiles  are  shown  completely  in  Appendix  III. )  A 
typical  velocity  profile  is  shown  in  Fig.  18. 

Numerically  integrating  such  functions  yielded  the 
coefficients  A^  and  Bn*  The  quality  of  the  curve  fit  is  determined 
by  how  well  it  approximates  the  original  function.  This,  in  turn, 
determines  the  number  of  coefficients  necessary.  In  obtaining 
Figs.  2,  22,  and  23  ten  coefficients  were  chosen  with 

-  mod  Bn)  =  VAn  - 

As  mentioned  earlier,  two  hot-wire  probes  were  used  so  that 
two  signals  could  he  recorded  simultaneously.  This  was  done  to 
obtain  data  for  the  auto-  and  cross-correlation  functions  associated 
with  the  flow  field* 

The  stationary  probe  was  used  as  the  reference  and  was 
located  between  two  fins  at  an  angular  position  ©  =  b^>  deg*  The 
axial  and  radial  position  of  both  probes  remained  constant  for  this 
test  at  0.95L  and  0-17D  respectively.  With  the  stationary  probe  in 
its  position,  the  movable  probe  was  traversed  around  the  body.  At 
selected  locations  of  A  ©  ,  the  angle  between  the  two  probes, 


the  signals  from  the  wires  were  recorded  on  the  magnetic  tape 
recorder  (CEC,  Type  GR  2800). 

This  method  was  used  to  obtain  data  for  both  the  b are-body 
configuration  and  for  the  body  equipped  with  propeller. 

Time  lags  were  introduced  by  the  computer  by  displacing  the 
data  during  the  computing  process-  To  be  able  to  use  the  available 
computer  facilities  for  the  actual  correlation,  the  data  had  to  be 
converted  to  digital  form.  This  meant  that  each  output  signal  (in 
volts)  was  transformed  into  an  integral  number  corresponding  to  its 
magnitude.  The  correlations  were  then  performed  with  these  numbers. 
Typical  correlation  functions  are  shown  in  Figs-  10  and  11. 

Mathematically,  the  auto- correlation  function  is  defined  by 

,T 

Ux«r)  U.  (ir  +  Cj  cbt  } 

-  T 

where  U is  the  signal  from  either  probe  representing  the 
velocity  fluctuations  measured  by  that  probe.  Note  that,  in  general, 
the  correlation  may  be  represented  by 
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where  x-j_,  x>>,  and  x^  are  the  displacements  along  the  coordinate 
axes-  Here,  only  is  considered-  (x^  denotes  the 

circumferential  direction,  y^  the  radial,  and  x^_  the  axial). 

Thus,  from  the  above,  the  autocorrelation  function  can  be 
obtained  for  each  probe  at  any  location. 
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In  general,  the  cross-correlation  function  is  defined  by 
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where  Uj  CJz)  and  U.^  (jt)  denote  the  fluctuating  velocities 
measured  by  the  probes* 

To  normalize  the  correlation  functions,  let  the  auto¬ 
correlation  for  zero  time  lag  of  the  movable  probe  be  denoted  by 
P\  and  that  for  the  reference  probe  by  •  Then, 

since  the  function  is  even,  the  limits  of  integration  can  be  changed 
to  0  and  T  with  the  value  of  the  integral  doubled.  Therefore, 
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Further  denote  the  autocorrelation  of  the  movable  probe  as  R 
and  that  of  the  stationary  probe  as  Ryy 
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For  the  cross- correlation  of  the  movable  probe  relative  to 
the  reference  probe,  denote  the  functions  as  for  and 

Ryx  ^or  ^  ^  O  - 

Using  the  previous  notation  for  the  autocorrelation  at 
r=o,  the  normalized  functions  are  given  by 
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The  cross-correlation  functions  in  their  normalized  form  are  shown 
in  Appendix  III  for  various  values  of  x*>. 

E»  Description  of  Procedure  and  Equipment  for  Obtaining  Correlation 
Functions  by  Computer  Methods 

The  computer  utilized  was  the  DDF- 19  made  by  Computer  Control 
Corp.  It  is  a  solid-state*  binary*  parallel*  19-bit-word* 
single-address*  general-purpose  computer.  The  memory  is  random 
access  with  a  magnetic  core  storage  of  4096  words.  The  machine 
access  time  is  3  microseconds  and  it  is  capable  of  100*000  additions 
per  second.  The  input-output  features  include: 

60-eharacter-per- second  paper-tape  punch 
300-character-per-second  photoelectric  paper-tape  reader 
15- character-per- second  typewriter 
300-line-pe.r-minute  printer. 

The  system  contains  two  magnetic-tape  recorders:  one  is  a 
direct-recording*  10- channel*  20-ineh-per- second  unit  with  a  packing 
density  of  333  bits  per  inch;  the  other  is  a  75-inch-per- second 
unit  with  packing  densities  of  200  and  556  hits  per  inch.  This  tape 
system  is  compatible  with  the  IBM  729  Mod  II  magnetic-tape  unit. 

Also  included  in  the  input-output  features  are  a 
200-increment-per- second  plotter*  three  digital- to- analog  conversion 
channels*  and  a  fast  serial  input  line  for  in- lab  testing  of 
experimental  systems. 

This  computer  has  one  built-in  index  register  of  12  bits  and 
an  interrupt  feature  for  priority  control. 


The  original  signal  consisted  of  the  hot-wire  output  and  was 


recorded  on  magnetic  tape  in  analog  form-  These  recordings 
represent  velocity  fluctuations  as  seen  by  the  wire  and  were 
proportional  to  some  voltage- 

The  tapes  were  then  processed  to  obtain  a  frequency  analysis 
of  the  signal-  These  analyses  give  an  indication  of  the  size  and 
concentration  of  the  energy  content.  A  typical  analysis  can  be  seen 
in  Fig.  2k -  Figure  25  shows  the  same  spectrum  but  from  wide-band 
frequency  considerations-  Once  these  tapes  were  so  analyzed,  it  was 
necessary  to  convert  the  signal  to  digital  form.  Because  of  the 
digitizing  equipment  involved,  the  original  tapes  had  to  toe 
re-recorded  to  filter  out  any  signal  above  1  kc-  This  was  done  by 
using  two  sets  of  Spencer -Kennedy  Laboratories  Inc-,  variable- 
frequency  electronic  filters,  Model  No-  J02,  and  resulted  in  the 
signal  shown  in  Fig  26-  This  figure  is  a  wide-band  analysis 
similar  to  Fig.  25-  The  rejection  rate  was  found  to  be  of  the  order 
of  50dB  per  octave  from  1  kc  to  2-kc- 

low  that  the  signals  were  revised  a  bias  voltage  was  applied 
in  order  to  bring  the  level  of  the  signal  to  a  range  from  0  to  10 
volts  d-c.  This  alteration  was  again  necessitated  by  the  electronic 
equipment  used  for  digitizing. 

With  the  voltage  now  in  this  range,  the  signal  was  digitized 
so  that  each  sample  voltage  corresponded  to  an  integer  between  1  and 
127«  This  is  the  range  used  by  the  computer  for  the  calculations  by 
binary  methods-  The  digitized  input  signal  was  sampled  at  a  rate  of 
3-3  kc-  This  corresponds  to  a  piece  of  data  1500  words  long.  (One 
word  is  a  seven-digit  number  representing  voltage  as  described 
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FREQUENCY  ~  (  kcps  ) 

Fig.  25  Wide-Band  Frequency  Analysis,  Before  Filtering,  for 
Same  Conditions  as  in  Fig.  24 


FREQUENCY  ~(  kcps  ) 

Fig.  26  Wide-Band  Frequency  Analysis,  After  Filtering  above  1  kc, 
for  Same  Conditions  as  in  Fig.  24 
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above. )  In  the  actual  sampling,  there  is  a  delay  of  300 
microseconds  between  words  so  that  the  sample  taken  amounted  to 
0.4-50  sec  of  actual  data.  It  was  subsequently  found  that  this 
sample  size  was  insufficient  to  give  repeatable  results.  A  method 
was  used  whereby  the  sampled  data  represented  2- 5  sec  of  actual  data, 
which  was  sufficient  to  give  repeatable  results. 

The  data  from  the  wires  was  recorded  on  two  channels  of  the 
CEC-GR2800,  seven-channel  tape  recorder. 

Figure  2J  is  a  sample  of  the  data  in  digitized  form-  Each 
jump  on  the  curve  represents  one  word.  When  calculating 
autocorrelation  functions,  the  signals  are  moved  relative  to 
themselves  and  relative  to  each  other  in  the  case  of  cross- 
correlation^  function  calculations.  Schematically,  the  relative 
displacements  can  be  shown  as  follows: 


I1500  words  -  wire  1 


t  >*- 


K - >|  1500  words  -  wire  lj - > 


for  the  autocorrelation,  and 


1500  words  -  wire  lj 


t 


- K - H  1500  words  -  wire  2 


for  the  cross-correlation. 


+  t 


+  t 


K  is  the  increment  by  which  the  data  are  displaced  by  the 
computer  when  performing  the  correlation-function  calculations.  It 
is  the  delay  time  of  300  microseconds  between  sample  words,  and  is 
an  integer  ranging  from  zero  to  $=1^^=1500.  As  seen  in  Fig.  10, 
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the  correlation  function  approaches  zero  as  K  approaches  100,  which 
corresponds  to  a  time  over  which  the  signals  are  correlated  of 
approximately  30  milliseconds* 

The  preceding  equations  for  the  required  correlations  did  not 
include  the  bias  necessary  to  raise  the  signal  level j  therefore,  it 
■was  removed  by  subtraction  when  the  correlation  functions  were 
computed.  This  was  accomplished  by  averaging  the  sample  words  and 
subtracting  this  average  from  each  word*  Once  the  words  were  so 
reduced,  the  actual  calculations  were  performed. 

For  any  given  time  delay,  K  =  0,1, ...I,  certain  words  from 
one  channel  line  up  with  words  from  the  other  channel  (see 
schematic).  These  combinations  were  then  multiplied  and  the  result 
divided  by  the  number  of  words  used,  I-K- 
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CHAPTER  IV 

THEORETICAL  CONSIDERATIONS  AND  DISCUSSION 

A.  Character  of  the  Solution 

The  theoretical  approach  to  the  problem  is  based  on  the 
solution  to  the  boundary- layer  equations  for  an  axisymmetric  body  of 
revolution  as  given  by  Truckenbrodt . (4) 

The  character  of  the  solution  is  such  that  only  two 
quantities  need  be  specified;  the  ideal  potential  flow  over  the 
original  body  and  the  Reynolds  number  based  on  body  length 

D 

With  only  these  quantities*  the  flow  field  can  be  obtained  for  the 
body  under  consideration.  This  flow  field  will  be  composed  of  a 
laminar  portion  to  transition*  transition  to  turbulent  flow  and* 
finally*  the  turbulent  region. 

The  solution  for  the  laminar  portion  up  to  the  transition 

/  j 

point  is  based  on  the  method  due  to  Von  Harman  and  Pohlhausen.  This 
solution  applies  to  two-dimensional  boundary  layers  and  is  used  with 
Truckenbrodt ' s  method  in  the  laminar  region.  According  to 
Truckenbrodt,  the  necessary  quantities  preceding  transition  are 
based  on  the  solution  for  the  two-dimensional  flat  plate.  His 
solution  applies  these  values  to  the  turbulent  boundary  layer  on  the 
axisymmetric  body  of  revolution.  This  completes  the  solution  for 
the  body  without  fins- 

The  fins  were  added*  and  the  solution  was  modified  to  account 
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for  their  presence.  Hie  effects  of  the  fins  on  the  flow  were 
determined  with  the  aid  of  Reference  7,  and  the  resulting  flow  into 
the  propeller  was  found- 

B-  Truekenbrodt 1 s  Method 

Truekenbrodt ' s  method  consists  in  using  the  energy  integral 
equation  to  calculate  the  momentum  thickness  rather  than  the  more 
familiar  momentum  integral-  The  basic  assumption  in  this  theory  is 
that  the  velocity  profile  in  the  boundary  layer  obeys  a  power  law  of 
the  form 

_u_  -  IJL\ /n 

u  'Is  J  '  a) 

The  basic  equation  for  the  distribution  of  the  momentum  thickness 
along  the  body  as  a  result  of  the  energy  integral  equation, 
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where 
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i/z  ^n+i)/n 
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is  a  constant  that  takes  into  account  the  laminar  portion  of  the 
■boundary  layer  up  to  the  transition  point  Xjt  / £  •  Here,  J[  is  the 

factor  by  which  linear  quantities  are  made  dimensionless. 

The  symbol  Cj.  in  equation  2  denotes  the  coefficient  of  ski: 
friction  of  a  flat  plate  at  zero  angle  of  attack  in  turbulent  flow 
at  a  Reynolds  number  =  — —  —  from  transition  to  X/jt  • 
The  symbol  in  equation  3  is  the  coefficient  of  skin  friction 

Ur,  X 


of  a  flat  plate  in  laminar  flow  at  a  Reynolds  number 


D 


from 


the  nose  to  the  transition  point 

The  quantity  01  j that  appears  in  both  equations  is  the 
ideal  potential,  velocity  distribution  for  the  body.  (This  is 
assumed  known  from  Fig.  28).  In  equation  2,  X  is  actually  a 
measure  along  the  body  curvature  rather  than  along  the  axis.  For 
convenience,  the  measure  of  curvature  has  been  changed  to  5  ; 
represents  the  axial  position.  Also,  R  j JL  is  the  radius  of  a 
cross- section  taken  normal  to  the  axis  at  the  position  X  jjl  « 

In  the  above  equations,  the  constant  D  has  been  empirically 
determined  and  is  given  as  6  for  turbulent  flow  at  a  large  Reynolds 
number.  In  this  case  J)  -  5*5x10°,  based  on  body  length. 

The  remaining  quantity  to  be  determined  is  the  variation  of 
the  shape  factor  L  is)  or  H ( S)  ,  This  will  then  give  the 
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required  displacement  thickness  §*  *  since  •  The 

shape  factor  (_  is  determined  by 


and  H  is  given  by  Fig.  29  once  the  value  for  L,  is  known. 

In  equation  4,  the  subscript  "  j  ”  refers  to  values  at  the 
transition  point.  The  symbol  is  essentially  a  variable  of 
integration  that  is  given  by 


where  H  =  6  and  fn  is  assumed  to  be  4.  NOTE:  at  transition* 


Once  the  values  at  transition  have  been  determined*  as  shown 
below*  the  momentum  thickness  is  found  from  equations  2  and  J. 

The  Reynolds  number  based  on  momentum  thickness  can  now  be 

found; 


57 


so  that  the 


quantity  bCf)  can  he  obtained  for  equation  4  as 


b(s)=o,o7  I°5ioR^-  0.23 


Thus,  LCS)  is  known  and  H  ( S )  is  found;  therefore;,  3*(S) 
is  specified. 

Co  Determination  of  the  Transition  Point-Method  of  Yon  Karman  and 
Pohlhausen  for  the  Laminar-Flow  Region 

The  problem  remains  to  find  those  quantities  associated  with 
the  laminar  flow  regime  and  the  determination  of  the  transition 
point  Xi/i  or  S*  •  The  original  body  configuration  is 
established  which  results  in  the  parameters  yj  and  yo/7  • 
The  potential  flow  solution  is  found  using  A-  M°  0.  Smith’s  method 
for  the  exact  solution  of  the  leumann  Problem.  (8)  Briefly,  the 
method  involves  the  numerical  solution  of  a  surface  source 
distribution  that  is  assumed  constant  over  finite  segments  of  the 
body-  The  coordinates  are  changed  from  axial  to  meridional  by 
simple  triangulation.  Thus, 


s 

Now,  knowing  the  relationship  between  U  j  Uao  and  ^5  >  the  Reynolds 
number,  based  on  momentum  thickness,  can  be  found  as  follows; 


(&jf U 


(8) 
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S  <  Si- 
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Following  PohUiausen ' s  method*  use  is  made  of  several  shape  factors 
such  as 


if  dU 

D  d  X 


(10) 


r  du 


^  d  X 


(ii) 


and  the  already  familiar  shape  factor*  f-j  =■  §  jD  »  In  the 
above*  S  is  the  actual  boundary- layer  thickness»  It  will  be  seen 
how  each  of  these  shape  factors  aids  in  finding  the  transition 
point • 


The  shape  factors*  K  and,_/\*  can  be  more  conveniently 
written  with  the  aid  of  equation  9  and  the  body  Reynolds  number . 
Thus* 


(12) 


It  has  been  found  that  and  _/\  satisfy  the  universal  relation 


K=  (a-i-bA  *cA2)a  , 


(13) 


where  CL  ,  b  ,  and  C  are  constants  whose  values  are  given  in 
Appendix  II.  The  shape  factor  H  can  also  be  written  in  terms  of 
K  and  .A; 
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where,  again,  cl  and  Q,  are  constants  and  are  also  given  in 
Appendix  II.  The  term  AH  is  an  empirical  factor  relating  the 
change  in  f-f  across  the  transition  region,  since  K  an<l  A.  aPply  only 
for  the  laminar  portion-  Here,  A,  H  was  chosen  of  the  order  of  1.2. 

The  necessity  for  the  above  development  is  to  find  the 
relationship  of  _/\.  and  S  in  terms  of  S  ,  the  arc  length  measured 
from  the  forward  stagnation  point.  Associated  with  the  flow 
downstream,  from  the  forward  stagnation  point  is  the  ability  of  the 
boundary  layer  to  remain  stable  and  not  to  become  turbulent.  The 
limit  of  stability  is  the  point  at  which  the  boundary  layer  changes 
from  the  stable  laminar  one  to  the  unstable  turbulant  one-  The 
Reynolds  number  associated  with  the  limit  of  stability  is  known  as 
the  "critical"  Reynolds  number.  As  long  as  the  local  Reynolds  number 
Urn  S*/l)  remains  less  than  the  critical  one  (u  rr)  $*Mcri+  ’ 
the  boundary  layer  remains  stable. 
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As  the  flow  proceeds  downstream,  the  rate  of  pressure 
decrease  becomes  smaller  and,  consequently,  is  followed  by  an 
increase  in  pressure  that  is  due  to  the  adverse  pressure  gradients 
Thus,  at  some  downstream  position,  Urp  §  /j)  increases  to  the 
local  limit  of  stability  (UmS'/V  )c  r  j  _j.  j  and  the  boundary 

layer  becomes  unstable.  The  position  associated  with  this  change  is 
the  transition  point;  it  is  defined  as 

s-stj  D  "l  S  Lit 


Certain  empirical  findings  permit  the  determination  of  the 
transition  point  with  relative  ease. 

Figure  JO  shows  the  variation  of  A  with  /umS*/i>lr;+  , 
which  then  gives  the  variation  of  £  with  (um  S*/-d)  i+  , 
since  £  vs  A  is  known  from  equation  11.  From  equation  9> 


vs 


S  is  known,  which  is  just  the  value  of  the  local  Reynolds  number 


since 
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Thus,  Fig.  J1  is  obtained,  which  shows  V ° /V){w  ^  l^jcr ii 
vs  S  "  By  the  definition  given  in  equation  15,  the  transition  point 
is  located  at  the  intersection  of  these  two  curves.  The 
displacement  thickness  S  and  momentum  thickness  A 
functions  of  S>  for  the  laminar  region  simply  as 


are  found  as 


Fig.  30 
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Fig.  31  Local  Reynolds  Number  and  Local  Limit  of  Stability  as  a 
Function  of  Position  along  Body 
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Therefore,  all  necessary  items  pertaining  to  the  laminar  flow  field, 
up  to  and  including  the  transition  point,  are  now  known.  ‘Thus,  the 
calculations  for  the  turbulent  field  may  now  be  carried  out  with  the 
above  knowledge.  The  values  at  transition  necessary  for  the 
remaining  calculations  are  given  in  Appendix  II. 

The  actual  location  of  the  transition  point,  as  seen  in 
Fig.  31,  agrees  with  data  presented  in  Reference  k  for  work  done  on 
an  8;1  ellipsoid  of  revolution. 

The  process  of  finding  the  displacement  thickness 
approximation  to  the  body  involves  the  use  of  equations  £  through  8 
in  a  process  of  iteration.  This  procedure  was  carried  out  until  the 

O  it  ,  -gj-j 

values  for  O  on  the  (n+l) '  approximation  were  within  5  per  cent 
of  those  found  for  the  n  "  case.  The  values  used  for  the  laminar 
portion  and  the  transition  region  were  those  found  after  the  first 
approximation  since  changes  in  these  quantities  are  negligible.  At 

the  rear  it  took  four  approximations  (ml=5)  to  satisfy  the  above 

0 

criterion.  The  final  values  chosen  for  a  ,  and  the  other 
quantities  associated  with  the  boundary  layer,  are  compared  with  the 
original,  body  coordinates  in  Appendix  II. 
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E-  Addition  of  Fins  to  the  Body 

After  the  solution  for  the  hare-body  configuration  had  been 
obtained,  the  fins  were  added.  Reference  7  is  an  experimental 
compilation  of  data  for  the  wake  profiles  of  struts  in  a  streaming 
flow-  With  the  inflow  velocity  at  the  strut  leading  edge  known,  the 
characteristic  parameters  associated  with  the  strut  wake  can  he  found. 

The  control- surface  profiles  are  MCA  0015  symmetrical 
airfoils  whose  profile  drag  coefficient  C<j0  was  obtained  with  the 
aid  of  Reference  9°  With  this  knowledge,  the  wake  depth  and  breadth 
of  influence  was  determined  (Fig-  32) » 

The  inflow  velocity  to  a  fin  was  found  from  the  boundary-layer 
profile  at  the  leading  edge  (Fig.  2-5,  Appendix  Il)»  This  velocity 
represents  the  actual  free- stream  velocity  seen  hy  the  fin.  The 
wake  profiles  were  obtained  at  the  corresponding  radial  distances 
used  in  the  experimental  measurements.  Thus,  at  a  particular  value 
of  r / D,  the  profile  is  determined  as  a  function  of  the 
circumferential  distance  6)  »  A  plot  of  the  profiles  thus  obtained 
are  shown  in  Fig-  3-  This  figure  compares  these  profiles  with  the 
measured,  ones- 

F.  The  Vortex  Effect 

Figure  3  shows  the  close  agreement  between  the  predicted 
wakes  and  the  measured  wakes  far  from  the  body- -such  as  at 
r/D  =  0.55“  However*  at  r/D  =  0.10,  which  is  close  to  the 
fin-and-body  intersection,  the  agreement  is  less  than  anticipated 
because  of  the  presence  of  corner  vortices  shed  by  the  fin  near  its 
root.  This  vortex  effect  is  the  only  flow  condition  that  had  not 
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been  accounted,  for  in  the  use  of  the  boundary- layer  equations.. 

Subsequent  testing  has  shown  the  presence  of  the  vorticity 
and  its  relative  strength.  Figure  4  shows  a  vortex  located  on  each 
side  of  a  fin.  The  vortex  was  measured  with  the  aid  of  a  vorticity 
meter  (Fig.  6)  developed  by  Dr»  !•  W.  McCormicks  Jr.  of  the 
Department  of  Aeronautical,  Engineerings  The  Pennsylvania  State 
University.  This  device  was  employed  as  a  means  of  establishing  the 
presence  of  the  streamwise  component  of  vorticity.  Figure  k  shows 
the  relative  strength  of  the  vortex  before  any  attempts  were  made  to 
eliminate  it.  If  the  vortex  could  be  eliminated,  the  propeller 
would  then  be  subjected  to  a  uniform  distribution  of  mean  velocity 
at  any  value  of  r/D  and  for  all  values  of  Q  (Fig.  7)°  In  this 
situation,  the  effects  of  the  unsteady  forces  would  be  favorably 
reduced. 

Efforts  made  to  reduce  the  strength  of  the  vortex  included 
a  technique  whereby  compressed  air  was  forced  out  of  holes  made  at 
the  fin-and-body  intersection  (Fig.  5).  Figure  7  shows  the  effects 
on  the  wake  resulting  from  an  air  pressure  of  5  psi.»  This  method, 
then,  appears  to  be  a  possible  means  of  eliminating  from  the  wake 
certain  characteristics  that  can  result  in  the  undesirable  operation 
of  a  propuls or. 


CHAPTER  V 


CONCLUSIONS  BASED  ON  EXPERIMENTAL  STUDIES 

A  Effects  of  the  Propeller 

The  results  of  the  circumferential  wake  distributions*  shown 
graphically  in  Appendix  III* are  summarized  in  Table  II. 


TABLE  II 

COMPARISONS  OF  RESULTS  OF  MEAN  AND  RMS  VELOCITY  BETWEEN 
MODEL  CONFIGURATIONS  AT  VARIOUS  VALUES  OF  r/D . 


r/D  Mean  Velocity  Ratio,  u/Uqo  RMS  Velocity  Ratio*;  V/Ucq 

Bare  Body  With  Prop  Bare  Body  With  Prop 


0.35 

0-9,33 

O.919 

0-0229 

0.0221 

0  •  24 

0.765 

O.876 

0-0483 

0.0241 

t— 

H 

6 

0-537 

0-697 

O.0691 

0.0447 

0.10 

0.326 

O.575 

0.0753 

0.0358 

(Note: 

the  design 

advan.ce  ratio  J  is  0-4o6 

for  the  model 

with 

propeller j  also,  r 

<  R  when  r/D  ^0-24) 

P 

A  comparison  of  the  measurements  shows  that  the  mean 


velocities  are  higher  for  the  body  with  propeller  than  without. 
This  is  true  when  the  radial  position  r/D  is  within  the  propeller 
radius  (r<Rp).  The  opposite  effect  is  noticed  when  traversing 
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radially  a  distance  greater  than  the  propeller  radius  (r>Rp).  It 
is  also  seen  that  the  rms  values  are  lower  with  the  propulsor,  which 
indicates  that  the  propulsor  has  a  smoothing  effect  on  the  flow* 
These  rms  results  have  agreed  with  previous  work,  and,  the  turbulence 
levels  obtained  peak  to  a  maximum  of  10  per  cent,  which  was  measured 
on  the  hare  body  for  r/D  =  0.10  (Fig-  3~9)«  It  is  also  to  be  noted 
that  the  mean-velocity  distribution  increases  in  magnitude  when 
traversing  radially  outward  from  the  body.  This  is  expected  because 
the  flow  field  is  greatly  influenced  by  the  turbulent  boundary 
layer,  which  behaves  in  the  same  manner-  This  influence,  which 
varies  inversely  from  the  body  outward  into  the  free  stream,  is 
negligible  once  the  flow  becomes  potential. 

It  is  felt  that  the  higher  mean  velocities  found,  with  the 
propulsor  are  due  to  the  streamlines  being  drawn  into  the  'body.  The 
propulsor  produces  this  effect  because  it  acts  as  a  favorable 
pressure  gradient  on  the  flow.  This,  then,  reduces  the  turbulence 
level  and  speeds  up  the  flow.  Hence,  the  thickness  of  the  boundary 
layer  is  decreased,  and  the  region  of  potential  f,low  is  greater. 

These  effects  produced  by  the  propeller  are  more  readily 
noted  in  the  boundary- layer  profile  than  in  the  wake  profile. 

Figure  33  clearly  illustrates  the  change  in  the  flow  as  it  passes 
through  the  propulsor.  The  change  in  the  flow  is  thought  of  as 
being  the  decrease  of  the  momentum  deficit,  which  results  from  the 
reaction  of  the  turbulent  boundary  layer  to  the  propeller. 

Figure  jk  shows  the  change  in  the  momentum  deficit  at  the 
propeller  inflow  region.  The  change  here  is  less  than  that 
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Fig*  35  Root-Mean-Square  Velocity  Distribution  for  Boundary-Layer 
Profile  at  Propeller  Inflow 


experienced  by  the  flow  as  it  passes  through  the  propuisor. 

Figure  35  further  illustrates  the  decreasing  turbulence  level  as  a 
result  of  the  smoother  flow  created  by  the  propeller..  This  figure 
consists  of  the  rms  values  for  the  boundary- layer  profile  shown  in 
Fig.  3h„ 

B-  Effects  of  the  Fins 

The  normalized  power  spectra,  obtained  by  passing  the 
hot-wire  signal  through  a  wave  analyzer  and  integrator,  are  shown  In 
Fig.  9-  These  spectra  are  for  the  model  with  a  propeller  at  two 
different  locations  of  r/D  and  &  •  From  Fig.  9  it  is  seen  that,  in 
each  case  chosen,  the  maximum  amount  of  energy  occurs  in  a  range  of 
nondimensional  frequency  ,  which  is  the  equivalent  of 

70-80  epg°  The  largest  amount  of  those  shown  occurs  at  the  closest 
position  to  the  body  and  midway  between  the  fins- 

Associated  with  this  energy  are  certain  eddies  produced  by 
the  shed  boundary  layer..  As  already  shown  by  the  Fourier  Analysis, 
the  effects  on  the  propeller  are  diminished  as  the  radial  distance 
increases  from  the  body.  By  noting  the  differences  between  Fig.  22 
and  Fig.  23,  and  also  from  Table  III,  it  is  apparent  that,  as  the 
radial  distance  traversed  exceeds  the  propeller  radius,  the 
influence  of  the  fin  at  the  inflow  region  becomes  negligible.  This 
was  noted  earlier  when  considering  the  region  of  potential  flow 
outside  the  influence  of  the  boundary  layer.  However,  it  is  seen 
that  the  flow  is  greatly  influenced  by  the  fins  when  traversing 
within  a  propeller  radius  of  the  body.  The  presence  of  the  large 
fourth  harmonic  close  to  the  body  (Fig.  22,  r/D  0.10  and  O0I7) 
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*Hotes  Under  each  value  of  r/D,  the  values  in  the  left-hand  column  were  obtained  with  all  fins  in  place; 
the  values  in  the  right-hand  column  were  obtained  with  only  the  forward  supporting  strut  in  place- 


gives  an  indication  of  the  disturbance  felt  by  the  blade  along  its 
span.  Figure  23  shows  that  the  blade  would  feel  relatively  no 
disturbance  along  its  span  in  the  absence  of  fins.  With  this 
knowledge,  the  propeller  designer  is  able  to  account  for  the 
unsteadiness  and  nonuniformities  encountered  by  the  blades.  Thus, 
the  response  associated  with  the  propeller-and-fin  configuration  is 
subject  to  control. 

C.  Harmonic  Content  Close  to  Body  Surface 

It  is  of  special  interest  to  consider  the  harmonic  content  of 
the  flow  field  in  a  region  close  to  the  body  which  affects  the 
propeller.  Figure  2,  and  the  accompanying  table,  compare  this 
content  at  a  value  of  r/D  =0=10  for  the  cases  of  the  body  and  fins, 
with  and  without  the  propeller.  From  Table  IV  it  is  seen  that  the 
magnitude  of  the  dominant  fourth  harmonic  is  reduced-  This  occurs 
because  the  induction  effects  of  the  propeller  reduce  the  influence 
of  the  turbulent  boundary  layer  by  smoothing  the  flow- 
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TABLE  IV 

COMPARISON  OF  HARMONIC  CONTENT  FOR  THE  STEADY  MEAN- VELOCITY 
PROFILE  AT  r/D  =  0.10  FOR  THE  BODY  WITH  FINS 


n 


Magnitude  of  the 
Fourier  Coefficient  Modulus,  Cn 

Without  With 

Propeller  Propeller 


1 

0.01713 

0.02330 

2 

0.04712 

0.00794 

5 

0.02638 

0.02340 

4 

0.14969 

0.11000 

5 

0.02134 

0.02060 

6 

0-00846 

0.00900 

7 

0.01378 

0.00392 

8 

0 . 02340 

0.01290 

9 

0 • 02001 

0.01740 

10 

0.01468 

0.00818 
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APPENDIX  I 

A  METHOD  FOR  CONDUCTING  REFERENCE  CHECKS  ON  THE  ELECTRONIC 
EQUIPMENT  USED  WITH  THE  HOT-WIRE  ANEMOMETER 

In  theory,  a  calibration  curve  of  output  voltage  versus 
velocity  should  be  a  straight  line.  However,  after  a  period  of 
time,  due  to  the  characteristics  of  the  electronic  equipment,  this 
straight-line  calibration  curve  could  change.  The  following  method 
allowed  the  hot-wire  output  to  be  checked  relative  to  that  of  a 
pitot  probe. 

In  a  wind-tunnel  facility,  other  than  that  being  used  for 
testing,  a  cylinder  was  mounted  normal  to  the  flow  and  the  hot-wire 
and  pitot  probes  were  positioned  behind  it.  The  pitot  probe 
establishes  the  mean  values  of  the  velocity  fluctuations  at  various 
tunnel  speeds  by  producing  a  difference  in  height  between  two 
mercury  columns.  From  these  pressure  readings  a  curve  was  drawn 
to  obtain  mean  velocity  versus  tunnel  speed,  which  shows  the 
response  to  the  turbulence  generator  (cylinder).  The  hot-wire  probe 
was  placed  in  relatively  the  same  position,  and  its  response  was 
recorded.  Once  the  mean  values  of  the  velocity  were  established 
with  the  pitot  probe,  any  necessary  changes  were  made  on  the  basis 
of  this  comparison  to  insure  the  proper  response  of  the  hot-wire 
probe. 

Once  it  had  been  established  that  this  hot-wire  probe  and  its 


accompanying  electronic  equipment  were  functioning  properly,  the 
prote  was  then  switched  to  the  electronics  for  the  other  probe .  The 
response  of  the  reference  probe  in  the  other  channel  was  observed 
and  necessary  changes  made  on  this  channel  to  insure  its  correct 
operation. 

By  means  of  this  checking  procedure,  the  hot-wire  anemometer 
was  then  used  to  obtain  quantitative  results  of  mean  velocity 
fluctuations  and  root-mean- square  values  of  the  turbulence  level. 
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APPENDIX  II 

RESULTS  OF  CALCULATIONS  FROM  THE  BOUNDARY - LAYER  EQUATIONS 
OF  VON  KARMAN  AND  POHLEAUSEN ,  AND  TRUCKENBRODT 

By  the  method,  outlined  in  Chapter  IV,  the  necessary 
quantities  associated  with  the  laminar  portion  of  the  boundary  layer 
were  obtained*  On  the  basis  of  Fig.  31,  the  transition  point  was 
obtained  satisfying  equation  15-  (These  computations  were  made  for 
the  body  without  fins  or  propeller. )  Table  2~I  presents  the 
approximation  to  the  body  found  for  the  laminar  displacement 
thickness  and  associated  quantities-  These  computations  were 
obtained  with  the  aid  of  an  IBM  1620  digital  computer.  As  seen 
from  Fig-  3 '-4  the  transition  point  occurs  for  S  »  0.124.5,  which 
corresponds  to  X/jl  =  0.111. 

The  associated  quantities  necessary  for  the  remaining 
computations  are  now  found  as  functions  of  the  transition  point; 

R  Rent  ~  1290  (from  Fig.  3D 

Cj  =  0-0033184  (for  equation  2) 

C ^  =  0-00170  (for  equation  3) 

1.034  (from  Fig.  28) 

H]  =  H  "A  H  =  1-322  (from£*?9  ,  equation  l4,  and  Reference  4) 
L,  =  0.0805  (from  Fig.  29) 

C*  =  0-1485  x  10"A  (from  equation  3) 
i?  s  (C*)~  0-4863  x  10"25  (from  equation  5) 
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TABLE  2-1 

APPROXIMATION  TO  ORIGINAL  BODY  BY  LAMINAR  DISPLACEMENT 
THICKNESS  FROM  THE  NOSE  TO  TRANSITION 

(All  dimensions  axe  normalized  with  respect  to  the  body  length. ) 


X 

Yo 

yp 

S* 

H-sy^ 

0,00000 

0-00000 

0.00000 

0.00000 

0.000000 

0.000 

. 00645 

.00949 

•00955 

.00004 

.000019 

2.100 

.02086 

. 02524 

•02551 

. 00070 

. 000055 

2.000 

. 05674 

•05045 

. 05054 

.00011 

. 000055 

2.000 

•05265 

•  05550 

. 05564 

. 00014 

. 000069 

2.000 

.06852 

•05955 

.05952 

. 00017 

. 000085 

2.000 

. 08425 

•  04250 

.04270 

. 00020 

. 000100 

2.000 

•10014 

•  04525 

. 04547 

. 00022 

. 000125 

1.760 

0 11100 

. 04668 

. 04695 

. 00025 

.000195 

1.522 

Note :  The 

last  value 

in  each  column 

is  for  the 

condition 

of 

transition. 
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The  constants  associated  with  equations  13  and  14,  which  are 
necessary  for  the  calculation  of  the  shape  factors  K  and  Jx *  are 
given  as  (4):; 


37 

315  ' 


1 

9^5 


_1_ 
9072 


and 


3 

10 


e  S3 


120 


At  transition,  K  and  jY  were  found  to  be  0.016  and  1.2* 
respectively. 

Once  these  values  at  transition  had  been  established*  it  was 
a  relatively  straightforward  procedure  to  account  for  the  turbulent 
portion  of  the  body.  The  use  of  equations  2  through  8  (section  IV) 
gave  successive  approximations  for  the  displacement  thickness. 

Table  2-II  and  Figs-  2-la  and  2-lb  show  these  approximations 
for  the  displacement  thickness  and  associated  quantities-  As  before* 
all  quantities  have  been  nondimensionalized  as  those  in  Table  2-1* 
with  y  being  the  approximation*  n  ~  1,2, 3*4. 

C*  7j 

In  addition  to  the  theoretical  predictions  for  O  and  (/ 
directly*  these  quantities  were  also  found  on  the  basis  of 
boundary- layer  profiles  at  the  probe  location.  In  section  III  0 
some  characteristic  results  of  test  measurements  are  shown: 
boundary  layer  profiles  are  shown,  Fig.  20,  and  S  and  'd  are 
defined. 

At  the  probe  position  xjjl  «  0.95,  S*  and.  2?  were 
obtained  from  Fig.  3 4*  without  propeller,  by  a  simple  integration 


Momentum  Thickness  for  the  Original  Body 


of  this  boundary” layer  profile.  The  integration  resulted  in 

n* 

o  =  0.00989, 

=  0.00551, 

and  H  =  1.80. 

By  interpolating  the  data  from  Table  2-II,  the  predicted  values  were 
obtained ; 

n* 

O  a  0-00935, 

'Q  =  0.00505, 

and  H  «  I.85. 

By  employing  equation  1  (section  IV )  with  H— 2=/f.H  “  1  j  > 
a  boundary- layer  profile  is  obtained  that  compares  favorably  with 
the  measured  profile.  Figure  2-2  shows  this  comparison  between  the 
measured  and  predicted  boundary  layers  at  the  probe  position. 

Figure  2-3  and  2-4  further  illustrate  the  use  of  equation  1 
for  obtaining  predicted  boundary- layer  profiles.  These  figures  are 
for  the  flow  at  a  strut  leading  edge  and  at  the  location  of  the 
propeller. 

Note;  The  apparent  "bulge”  in  the  dotted  curve  of  Fig.  2-2 
shows  that  the  flow  is  starting  to  reverse  itself,  which  is  an 
indication  that  the  separation  point  is  at  hand.  This  effect  is 
also  noticeable  in  the  predicted  curve  of  Fig  2-2,  and  in  Fig.  2-4, 
which  appear  as  a  "shallowness",  when  compared  with  Fig.  2-3,  which 
is  a  "full"  boundary- layer  profile. 


APPROXIMATION  TO  ORIGINAL  BODY  BY  DXSPLA< 
THICKNESS  FROM  TRANSITION  TO  TAIL 
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Pig.  2-4  Predicted  Boundary- Layer  Profile  at  Propeller  Position 
(From  Equation  1) 


APPENDIX  III 


GRAPHICAL  PRESENTATION  OF  THE  EXPERIMENTAL 
MEASUREMENTS  MADE  IN  THE  BOUNDARY  LAYER 

The  complete  results  of  the  hot-wire  measurements  are 
presented,  in  Figs-  J-l  through  3-9-  These  results  include  the  mean 
velocity  and  rms  velocity  distribution  as  a  function  of 
circumferential  angle  of  traverse  9»  Each  graph  is  for  a  different 
value  of  .r/D  for  the  test  conditions  of  the  body  and  fins  both 
with  and  without  propeller.  Figure  3~1  gives  the  mean  velocity 
distribution  for  all  values  of  r/D  on  one  graph  to  show  the  overall 
variation  with  distance  from  the  body.  Table  II  was  established  on 
the  basis  of  the  average  values  of  these  results- 
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Fig-  3-1  Steady  Mean-Velocity  Distribution  as  a  Fun- 
Radial  Distance  for  Body  without  Propeller 


Fig-  3-3  Root-Mean-Square  of  the  Random-Velocity  Fluctuations 
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Fig.  3-6  Steady  Mean-Velocity  Distribution  for  r/D  =  0.17 


Fig.  3-7  Root-Mean-Square  of  the  Random- Velocity  Fluctuations 
for  r/D  =  0.17 


Fig.  3-8  Steady  Mean-Velocity  Distribution  for  r/D  =  0.10 


Fig*  5-9  Root-Mean-Square  of  the  Random-Velocity  Fluctuations 
for  r/D  =  0.10 
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The  following  figures  are  results  of  the  cross  correlations 
obtained  at  representative  distances  of  probe  separation..  These 
cases  are  for  the  body  with  fins.  In  Figs.  3-10  and  3-11  the 
distance  between  the  probes  is  0.25  in.  Figures  3-12  and  3-I3  are 
for  a  distance  between  probes  of  5-9  in.  These  figures  represent 
the  extreme  cases  of  probe  separation. 

In  addition  to  these  cross-correlation  functions,  their 
associated  cross- spectral,  densities  are  also  shown.  As  before, 
these  were  obtained  by  performing  a  Fourier  transform  of  the 
original  functions.  The  cross- spectral  densities  are  shown,  in 
Figs.  3-14  through.  3-21- 
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Fig.  J>-11  Cross-Correlation  Function  for  Body  and  Fins  without 
Propeller  for  =  0.25  inches 
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Fig.  3-12 


105 


X 

>> 


Ct 


>* 

X 

ce 


UJ 

o 

Li_ 

Ll 

LlI 

o 

o 


o 

I- 


Fig.  3-13 


Cross- Correlation  Function  for  Body  and  Fins  without 
Propeller  at  x^  =  5*9  inches 
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Fig.  3-15  Imaginary  Part  of  Cross-Power  Spectral  Density  for  Body 
and  Fins  with  Propeller  for  x,  =  0.25  inches 
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Fig.  3-16  Real  Part  of  Cross-Power  Spectral  Density  for  Body  and 
Fins  without  Propeller  for  x  =  0.25  inches 
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Fig.  3-17  Imaginary  Part  of  Cross-Power  Spectral  Density  for  Body 
and  Fins  without  Propeller  for  x  =  0.25  inches 
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Fig.  J-lQ  Real  Part  of  Cross-Power  Spectral  Density  for  Body  and 
Fins  with  Propeller  for  xx  =  5*9  inches 
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Pig.  3- 19  Imaginary  Part  of  Cross-Power  Spectral  Density  for  Body 
and  Fins  with  Propeller  for  x,  =  5*9  inches 
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Fig.  3-20  Real  Part  of  Cross-Power  Spectral  Density  for  Body  and 
Fins  without  Propeller  for  x  =  5*9  inches 
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Pig.  3-21  Imaginary  Part  of  Cross-Power  Spectral  Density  for  Body 
and  Fins  without  Propeller  for  x  =  5*9  inches 
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